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Abstract. The results of a new spectroscopic survey of 66 z £ 4 quasars for 
Damped Lyman-a absorption systems are presented. The search led to the discovery 
of 30 new DLA candidates which are analysed in order to compute the comoving 
mass density of neutral gas in a non-zero A Universe. The possible sources of 
uncertainty are discussed and the implications of our results for the theories of 
galaxy formation and evolution are emphasized. A subsequent paper will present 
details of the calculations summarised here and a more extensive explanation of the 
consequences of our observations for the understanding of the nature of DLAs. 



1. A New Survey for Damped Lyman-a Galaxies 

Damped Lyman-a (hereafter DLA) galaxies are defined as absorption 
systems in the spectra of quasars with a hydrogen column density 
N(HI) > 2 x 10 20 atoms cm -2 (> 10 A rest frame equivalent width, 
Wolfe et al. 1986). Their space density distribution, metal abundances 
and kinematics seem to indicate that they are the progenitors of present 
day galaxies, thus providing a powerful observational tool for galaxy 
evolution studies. 

DLAs probe the neutral gas from which stars form and thus are a di- 
rect means to trace galaxy formation. Quasars are extremely luminous 
objects which are observed up to z = 5.8 (Fan et al. 2000), and thus 
DLAs can be probed over a large redshift range. Finally DLA galaxies 
can be detected to fainter limits than objects detected using traditional 
emission observations. 

We observed 66 z ^ 4 quasars (see Peroux et al. 2000 for details on the 
observations) from the 2 nd APM colour survey (Storrie-Lombardi et al. 
2000), the Palomar Sky Survey (Kennefick et al. 1995a & 1995b and 
G. Djorgovski's web page) and the Sloan Digital Sky Survey (Fan et 
al. 1999). This resulted in 31 DLA candidates, 30 of which are newly 
discovered and many with associated metal lines. 
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Damped Lyman— a Absorption Redshift 

Figure 1. Sensitivity functions for DLA galaxy searches. The dotted line is from 
Storrie-Lombardi et al. 1996 (A(z) = 239), the dashed line is the compilation of 
Storrie-Lombardi & Wolfe 2000 (A(z) = 418), and our survey is represented by 
the solid line (A(z) = 493). Our new sample more than doubles the redshift path 
surveyed above z £ 3.5. 



2. Probing the Neutral Gas Content of the Universe 

For each quasar the redshift range along which a DLA could be observed 
is determined and used to deduce the survey's sensitivity. This is repre- 
sented by the g(z) function shown in Fig 1. Our new data is analysed in 
conjunction with previous observations resulting in 115 DLAs observed 
in 697 quasars (Rao and Turnshek 1999 data not included). 
The mass of neutral gas contained in DLAs is expressed as: Q,dla{z) = 
H °f™f f™ m . Nf(N, z)dN, where N is the column density and f(N, z) 
is the distribution function. The integral is estimated by summing up 
over all the observed column densities and dividing by the absorption 
distance interval, X(z), which corrects for comoving coordinates and 
thus depends on the geometry of the Universe. In a non-zero A Universe: 
X(z) = j z (l + z) 2 [(1 + z) 2 {l + zQ M ) ~ z(2 + z)n K y 1/2 dz. 

Our observations are presented in Fig 2 for a currently popular 
cosmology and show the evolution of Qdla over the last 90 % of the 
age of the Universe thus probing the early stages of galaxies even before 
gas turns into stars. 

In order to account for the neutral gas not included in DLAs, we used a 
power law with an exponential turnover (similar to the Schechter func- 
tion) in extrapolating from the observed column density distribution: 
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Figure 2. The circles show the neutral gas in Damped Lyman-a galaxies (the filled 
circles are our observations and the open circles are Rao & Turnshek 1999 results) . 
Vertical bars correspond to 1 a errors and the horizontal bars indicate bin sizes. The 
stars are the total HI including a correction for the neutral gas contained in systems 
with column densities below 2 x 10 2() atoms cm" 2 . The triangle at z = is from 
Natarajan & Pettini 1997 who used a recent galaxy luminosity function to reproduce 
the local HI mass measured by Rao & Turnshek 1993. The squares, Q,fhp and flao 
(Fukugita, Hogan & Peebles 1998 and Gnedin & Ostriker 1992 respectively) arc 
Qbaryons in local galaxies. The semi-analytical models from Somerville, Primack & 
Faber 2000 are overplotted (see text for detailed explanations of the models). 



f(N,z) = (f*/N*)(N/N*) @e N I N *. Using the parameters determined 
by Storrie-Lombardi et al. 1996, iV* = 21.63 ±0.35 and (3 = 1.48 ±0.30, 
results in a ~ 20% increase in ft for z < 3.5. For z > 3.5 systems, our 
observations indicate a steepening of the low column density end of the 
DLA distribution (where much of the mass resides) in comparison with 
our previous study and leads to a factor of ~ 2 correction (as shown 
by the stars in Fig 2). 



3. Discussion 



Most of the uncertainties discussed here seem to indicate that £Idla 
is being underestimated. At low redshift, Rao & Turnshek 1999 used a 



peroux-c.tex; 1/02/2008; 14:52; p. 3 



4 



Peroux et al. 



different method (based on the observations of Lyman-a in identified 
Mg II systems) and HST data to derive a higher neutral gas content 
in DLAs (Fig 2). These results are surprisingly high and might be due 
to low number statistics or might indicate that high redshift results 
are underestimated. More uncertainties arise in the fact that dust is 
likely to produce observational biases unaccounted for here and that 
the exact nature of DLAs is unknown. In particular, they might be 
different types of object at different redshift. Finally, higher resolution 
data are necessary in order to identify multiple systems. 
Several groups (Kauffmann & Haehnelt 2000, Somerville, Primack & 
Faber 2000) have included more realistic physics in their simulations 
to construct semi-analytical models of galaxy formation which, among 
other things, predict the evolution of cold gas in the Universe. The 
models presented in Fig 2 vary in their recipe for star formation: it is 
triggered by galaxy-galaxy mergers in the collisional starburst model, 
constant with redshift in the constant efficiency model, and scales in- 
versely with disc dynamical time in the accelerated efficiency model. 
Our observational results can thus be used to directly constrain theories 
of galaxy evolution. 

We would like to thank Isobel Hook for help in acquiring the data. 
CP is grateful to Max Pettini and the organisers for help with a very 
enjoyable conference. 
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